Influence of electron-acoustic phonon scattering on off-resonant cavity feeding within 

a strongly coupled quantum-dot cavity system 

S. HughetB and P. Yao 

Department of Physics, Queen's University 
Kingston, ON K7L 3N6 Canada 



O 

(N 

m 



O 

o 



> 

(N 
t> 
O 

(N 
O 



X 



F. Milde and A. Knorr 

Institut fur Theoretische Physik, Nichtlineare Optik und Quantenelektronik, 
Technische Universitat Berlin, Hardenbergstrafie 36, PN 7-1 10623 Berlin, Germany 

D. Dalacu, K. Mnaymnch, V. Sazonova, P. J. Poole, G. C. Aers, J. Lapointc, R. Cheriton, and R. L. Williams 

Institute for Micro structural Sciences, National Research Council, Ottawa, Canada K1A 0R6 

We present a medium-dependent quantum optics approach to describe the influence of electron- 
acoustic phonon coupling on the emission spectra of a strongly coupled quantum-dot cavity sys- 
tem. Using a canonical Hamiltonian for light quantization and a photon Green function formalism, 
phonons are included to all orders through the dot polarizability function obtained within the inde- 
pendent Boson model. We derive simple user-friendly analytical expressions for the linear quantum 
light spectrum, including the influence from both exciton and cavity-emission decay channels. In the 
regime of semiconductor cavity-QED, we study cavity emission for various exciton-cavity detunings 
and demonstrate rich spectral asymmetries as well as cavity-mode suppression and enhancement 
effects. Our technique is nonperturbative, and non-Markovian, and can be applied to study photon 
emission from a wide range of semiconductor quantum dot structures, including waveguides and 
coupled cavity arrays. We compare our theory directly to recent and apparently puzzling exper- 
imental data for a single site-controlled quantum dot in a photonic crystal cavity and show good 
agreement as a function of cavity-dot detuning and as a function of temperature. 
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I. INTRODUCTION 

The influence of electron-acoustic phonon scattering is 
a well known effect in semiconductor quantum dot (QDs). 
Due to the interplay of phonon emission and absorption, 
the longitudinal acoustic (LA) phonon bath manifests in 
spectral lineshapes that are highly asymmetric at low 
temperature, and sit on the background of the symmet- 
ric zero phonon line (ZPL). The characteristic spectral 
lineshape of LA phonon scattering on electron-hole pairs 
( "excitons" ) in quasi- homogeneous semiconductor struc- 
tures is well described through the independent Boson 
model (IBM)i~— , and IBM simulations have shown good 
agreement with experiments. In addition to the IBM line- 
shape, the broadening of the ZPL is usually described 
phenomcnologically, though this can be reliably fit to 
experiments^; however, there is still some controversy as 
to the origin of the ZPL broadening which can include 
contributions from radiative broadening, spectral diffu- 
sion, anharmonicity effects^, phonon scattering from 
interfaces^, and a modified phonon spectrum^. While 
it is well known that the phonons cause the exciton line- 
shapes to be highly non-Lorentzian (stemming from non- 
Markovian decay), most of the present QD cavity-QED 
theories only add in a Lorentzian broadening mechanism 
for the QD excitons, e.g., at the level of a Markovian 
master equation. 

Recently, there have been several studies of the role 
of electron-phonon coupling in scmicondutor QD cav- 



ity systems. Wilson-Rae and Imamogliiii treated the 
phonon interaction with QDs using polaron Green func- 
tions and derived an analytic absorption lineshape when 
the dot and cavity are on resonance; in the polaron rep- 
resentation, new phonon-induccd interaction terms are 
introduced exactly, while a second-order Born approxi- 
mation was applied to include "residual" exciton-photon- 
phonon coupling effects. Polaron approaches have also 
been recently employed to successfully describe phonon- 
induced decay of optical pulse-excited QDs in the ab- 
sence of any cavity coupling^. Milde et al*£ numerically 
solved the IBM and coupled the QD susceptibility to a 
photonic crystal cavity system through a scmiclassical 
Green function approach, demonstrating asymmetries in 
the on-resonance Rabi doublet and the effect of increasing 
temperature for both the cavity emitted spectra and the 
side-coupled waveguide transmission. Xue et aL— have 
applied perturbation theory to study the phonon-induccd 
decoherence on vacuum Rabi oscillations as a function of 
detuning between the cavity mode and exciton. Kaer et 
al*£- also explored off-resonant interactions between the 
QD and a cavity using a numerical solution to the sys- 
tem master equations within a time-convolutionless ap- 
proach. Ota et alr^- numerically solved the Wilson-Rae 
and Imamoglii system master equation^, and demon- 
strated the importance of asymmetric off-resonance cou- 
pling and non-Markovian relaxation, finding good agree- 
ment with their experiments. Hohcnester— introduced a 
useful model to derive an effective phonon-mediated scat- 
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tcring rate (see also Refill), and has found good agree- 
ment with experiments in the weak coupling regime^ 8 -. 

From an experimental viewpoint, significant off- 
resonant coupling between an exciton and a cavity has 
been seen in a number of semiconductor QD-cavity sys- 
tems, e.g., see RefsJ^ - — ; and several theoretical and ex- 
perimental works have tried to explain the basic cou- 
pling mechanisms^ - — , most of which stem from the ex- 
citon broadening mechanisms and the simple physics of 
two coupled oscillators. Given the importance of ex- 
citon decay processes on the QD-cavity coupling, one 
can therefore expect that the interaction of phonons can 
play a qualitatively important role when certain coupling 
conditions are met. However, while the aforementioned 
theoretical phonon studies are interesting in their own 
right, unfortunately none of them present simple analyt- 
ical spectra that allow one to explore a wide range of 
coupling phenomena both on- and off-resonance with a 
strongly coupled cavity. Thus many of the experimental 
- and theoretical - groups continue to use the only avail- 
able analytical formulas with no phonon coupling directly 
included, e.g., see Ref<22. 

Here wc apply a photon Green function method to de- 
rive useful analytical spectra, with phonons included to 
all orders through the exciton polarizability. Although a 
similar semi-classical Green function approach was pre- 
sented by Milde et alJ£, only on-resonance conditions 
were studied, and the IBM was solved numerically. Tarel 
and Savona3 have recently presented a semi-classical 
Green function spectra with phonons included to sec- 
ond order (also previously presented in Ref^), where 
numerical solutions of the phonon baths were exploited. 
We first study leaky cavity emission with and without 
phonons and connect the results to the recent work of Ota 
et al^, and find good qualitative agreement with their 
observations, namely pronounced asymmetries for high 
energy or low energy cavity coupling and an asymmet- 
ric vacuum Rabi doublet (at low temperatures, T » 4 K); 
we also demonstrate that phonon-induced cavity suppres- 
sion can occur, which is an effect that stems from the real 
part of the phonon self-energy. Secondly, we compare di- 
rectly with recent experimental measurements by Dalacu 
et al^, who studied QD-cavity coupling in single site- 
controlled QDs and found that the data could not be fit 
without adding in some unknown (detuning-dependent) 
cavity pump term; in contrast, we show that this data 
can be well reproduced using our analytic model with no 
cavity pump term included at all. Importantly, our spec- 
tra can be applied to a wide range of systems, including 
waveguides, and it rigorously applies to both weak and 
strong coupling regimes, and contains phonons to all or- 
ders (at the level of the IBM) . We also show higher tem- 
perature results and systematically compare these with 
experimental data. In our experimental-theory analysis, 
we include two decay channels, accounting for both the 
radiation-mode emission and the leaky cavity emission. 

Our paper is organized as follows. In Sec. UH we intro- 
duce the basic theory and analytical formulas for calcu- 



lating the emitted spectra from a strongly coupled QD- 
cavity system in the presence of electron-acoustic phonon 
coupling. Example computed emission spectra with and 
without phonon coupling are shown in Sec. IIIII Sec- 
tion IIVI compares our theory with recent data on site 
controlled single QDs in photonic crystal cavities^, and 
demonstrates the significant influence of electron-phonon 
coupling. In Sec. E] wc conclude. 



II. THEORY AND ANALYTICAL SPECTRA 

For this work, we are interested in the linear spectrum 
and thus consider a QD-cavity system that is weakly 
pumped, incoherently, where the emission dynamics stem 
from an excited electron-hole pair. In general, there have 
been several theoretical approaches to this problem in the 
literature. Two of the most powerful methods include 
Green function approaches^ 8 - and quantum master equa- 
tion techniques 3 -! - — . A major advantage of the Green 
function approach is that one can obtain analytical spec- 
tra for any inhomogeneous and lossy structures, including 
lossy metamaterial waveguides 3 -^ and a variety of coupled 
cavity- waveguide system a 35 ' 36 . In contrast, a key advan- 
tage of the master equation approach is the ease with 
which it adds in additional dissipation effects such as 
pure dephasing, though at the level of obtaining emis- 
sion spectra, both Green function and master equation 
approaches can be equivalent. A significant disadvantage 
of the master equation approach is that it is typically 
limited to simple leaky cavity systems with a Lorentzian 
decay, i.e., Markovian theory. Master equation solutions 
can also consider more realistic initial conditions, such as 
those obtained through steady-state pumping. 

The general Green function approach to obtain- 
ing the electric-field operator has been described 
elsewher o 28 ' 37 i 38 . Here we will briefly highlight the the- 
oretical background and concentrate on presenting the 
general expressions for the field operator and the spec- 
tra. At first we neglect non-radiative broadening on the 
exciton decay, which allows us to obtain the exact analyt- 
ical field operator without coupling to phonons. Specifi- 
cally, we use a canonical Hamiltonian that quantizes the 
macroscopic electromagnetic fields, and exploit the dipolc 
approximation for the QD-mcdium coupling: 

H = fiio x a x a x + huj\a\a\ 

A 

-ih£2(a- + &+)(g x a x - g* x a{), (1) 

A 

where a\ represents the field mode operators and 
are the Pauli operators of the QD excitons, lo\ is the 
cigcnfrcquency corresponding to the transverse modes of 
the system [f\(r)], excluding the dot; g\ is the field-dot 

coupling coefficient, defined through g\ = y^f^Mx • 

fx(r), with (j, x = n x [i x the optical dipole moment of 
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the exciton, aligned along n x (a unit vector). We con- 
sider only one target exciton in the spectral region of 
interest for the coupled QD and assume that coupling 
to the other polarized exciton is negligible. The Heisen- 
berg equations of motion for the operators can be used 
to derive the electric field operato r 28 ! 38 . Considering a 
weak excitation condition (i.e., we neglect higher-order 
photon-correlation effects, which is valid in these sys- 
tems for weak powers^ 3 -), and assuming an excited QD 
in vacuum, we derive the quantum field operator—: 



E(r,u;) 



G(r, r d ;uj) ■ d x (uj)/e 
1 - n x ■ G(r d ,r d ;uj) ■ n x a x (iv)'' 



(2) 



where is the QD position, a x {ui) = ^f"^ is the 
bare (no radiative or non-radiative coupling) exciton po- 
larizability, and d x (uj) = —ifi x [a~(t = 0)/(u) — u) x ) + 
(r x {t — 0)/(lj + u) x )] is a quantum dipole source that 
originates from the excited QD. Without phonon inter- 
actions, the above operator is exact within the stated 
model approximations. The propagator, G(r, r';w), is 
the transverse^ 3 , photon Green function of the medium, 
without any QD, and for cavity systems it can be easily 
written down analytically in terms of the cavity mode de- 
cay (cavity emission) and the radiation-mode decay (dot 
emission)^ 8 .. Although electric field pumping (e.g., cav- 
ity pumping) can be included in a straightforward way, 
we are primarily interested in the spectra from an ex- 
cited (and single) QD, so we treat the initial field as in 
vacuum. 

With a knowledge of the medium-dependent Green 
functions, and a suitable initial condition for exciting the 
material system, one can conveniently obtain the ana- 
lytic emitted spectrum, in any general structure, from^ 8 . 
S(r,w) = ([E(r, cj)]tE(r, cj)). The above expressions are 
exact, and no assumption has been made about the form 
of the medium. 

For a semiconductor cavity system such as a pla- 
nar photonic crystal cavity, one easily obtains the cav- 
ity emitted spectrum and the radiation-mode emitted 
spectrum 2 ^, analytically, 



S r (r,w) = F r (r)r rad x 

LO x + LJ 



S c (r,ij) =F c {r)T c x 

2gui c {uj x 
uj 2 — uj 2 —iujT c 



( ,2 , ,2 ,p 4:g 2 uj x uj c 

x x ujt—uj 2 —iujY c 



(3) 



(4) 



through Y x = r rac j + V , with T' due to pure dephasing 
processes; in general, for cw (continuous wave) pumping, 
both phonon effects and spectral diffusion will enhance 
r', and r' is also known to be temperature dependent^. 
It is important to note that in a planar photonic crystal 
structure, both S r and S c photon decay channels con- 
tribute to vertical photon emission. The radiation-mode 
decay is due to coupling to the continuum of radiation 
modes above the slab light line^; and the F c / r (r) rep- 
resent the geometrical factors which will depend upon 
the collection geometry of light emission. For a micropil- 
lar cavity system 3 ^, typically only the cavity emission is 
required in an identical form to above, and so the pre- 
scriptions apply to a wide range of semiconductor cavity 
systems. 

To include phonon interactions in a simple but rigorous 
way, we assume that the cavity and phonon correlation 
functions can be decoupled, and add in the phonon polar- 
izability via the known phonon self-energy, T, p h(to), from 
the IBM In essence, we are then considering the opti- 
cal polarizability of the QD in the presence of phonons, 
as the exact perturbation to the medium. Consequently, 
one has a slight modification to the above spectra, 



Sr(r,w) = F,.(r)r rad x 

UJ X +LU 



2 X -UJ 2 - iujT x - wE p h(w) - u -i-J-t,T e 
S c (t,u) = F c (r)T c x 



2guj c (lj x -\-uj) 

— 2 7 5 : — 



id 2 — IUJT T — UjYi n h(uj) 4 g ^x^, 



, (5) 



, (6) 



where the cavity emission is similar in form to the one 
presented by Tarel and Savona^, where phonons were 
included to second-order and a rotating-wave approxi- 
mation was made. As limits, we obtain the correct IBM 
spectral form for exciton decay, and earlier derived spec- 
tra for semiconductor cavities^. We also obtain ZPL 
broadening associated with the leaky cavity system. 

In the spirit of deriving a simple analytic solution, with 
phonons included to all orders, the strategy is to use an 
analytic phonon self-energy at the level of the IBM. To 
do this, we exploit phonon spectral functions, similar to 
the ones used by Wilson-Rac and Imamogliiii, but we 
use a more appropriate spectral function for phonon in- 
teractions via a deformation potential^; similar spectral 
function are commonly used when describing LA-phonon 
coupling. To obtain the phonon self-energy, the IBM 
time-dependent phase must be added into the Lorentzian 
decay model for the exciton; this is obtained from 



with the total spectrum S t (r,uj) = S r (r,u)) + S c (r,ui). 
We emphasize that radiative coupling to the cavity sys- 
tem is fully included by coupling to both the contin- 
uum of radiation modes and the leaky cavity mode — 
with a decay rate given by r c . Additional broadening 
of the ZPL has also been included, phenomcnologically, 



duiJ(u>)/u} 2 \coth.((3hu}/2) cos(wt) — isin(wi)] 



(7) 

that describes the electron-LA-phonon interaction, at 
temperature T = l/(3kb; we consider electron-phonon in- 
teractions via a deformation potential coupling which is 
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known to account for the major phonon interactions in 
our considered QD. Although one may have LO inter- 
actions as well, these give no dephasing. We consid- 
ering a spherical QD model, with similar electron lo- 
calization lengths in the valence and conduction bands 
(l e = 1^ pa 5nm). Generalizing to include electrons 
and holes with different localization lengths is straight- 
forward, but the equations become more cumbersome. 
Thus, spectral function can be conveniently defined as 4 ^ 



J(u>) = a p ui 3 exp(— uj 2 /2lo 2 ), 



(8) 



where we use Wf> = 1 meV and a p j (2tt) 2 = 0.06 ps 2 as rep- 
resentative numbers for InAs type QD o 14 ' 41 . The defor- 
mation coupling constant used here is somewhat smaller 
than the value used in Refj 1 ^ (where a p / (27r) 2 w 0.1 ps 2 ), 
though in the literature there is no well accepted values 
for InAs QDs; moreover, the dimensionless Huang-Rhys, 
Shr = a p /{2n) 2 'cf/lK « 0.034 (where c x = 3800m/s 
is the speed of sound), has been shown to be signifi- 
cantly enhanced in QDs (from its bulk value); for ex- 
ample, for InAs/GaAs QDs, S H r = 0.01 - 0.5 (Ref42) 
and Shr = 0.5 (Ref4£) have been reported, and vari- 
ous mechanisms for such enhancements have been pro- 
posed, including defects and non-adiabatic effects. The 
phonon interactions also result in a polaron shift, A = 
/ °° dui J {lu) / uj = Shr « 42 fxeV. In what follows 

below, we neglect the polaron shift as it merely adds a 
fixed frequency, and we can always redefine u> x , so our u> x 
includes the polaron shift. We also point out that there 
will likely be other non-diagonal phonon coupling as well 
and so in principle the phonon parameters above could 
be varied and used to fit experiments. 

With the above analytical form for the phonon bath, 
the time-dependent polarizability takes the form a x (t) = 
a x (0) exp[—i(w x + A)t + ip(t)], an d t ne frequency de- 
pendent polarizability is obtained from a simple Fourier 
transform, yielding 



a x (io) 



d 2 2uj x /heo 



iloT x - ujT, ph (uj) 



(9) 



To help better explain the phonon coupling effects shown 
later, in Fig. [T] we show the phonon self-energies for the 
two different temperatures of T = 4 K (a) and T = 40 K 
(b). We recognize a local minimum in Im[Ep^] (red-solid) 
near uj x , a maximum in Im[E p /j] near lmeV, and signifi- 
cant Re[E p h] (blue-dashed) over a broad spectral range. 
These results imply that both the real and imaginary 
contributions will have a significant impact on phonon 
coupling effects in the regime of cavity-QED. 

As mentioned above, our general methodology is sim- 
ilar in spirit to a semi-classical approach. Specifi- 
cally, we have assumed a well defined spectral line- 
shape for the QD susceptibility, and then coupled this 
QD frequency response - with phonons included self- 
consistcntly - to the medium-dependent Green func- 
tions to obtain the analytical spectra. Comparing 
with the approach of Wilson-Rac and Imamoglu, they 



have in their system Hamiltonian a phonon-modificd 
cavity coupling rate 1 - 1 - g — > g{B), with (B) = 
cxp(— 0.5 J dujJ{ui)coth(Phu}/2)/uj 2 ), which, using the 
parameters above for T = 4 — 40 K, is around (B) = 
0.91 — 0.55. We do not have this term explicitly, however 
our self-energy term naturally includes this coupling, and 
to all orders. To make this clearer, we can rewrite the 
solution as, e.g., for the cavity- mode emission, 



S c (r,w) =F c (r)r c x 



>-u' 2 -iu)r T 



ioST r 



(10) 



where the real part of the phonon self-energy causes a 
reduction in the g coupling. Such an effect was previously 
pointed out before for on-resonance coupling between a 
cavity and an exciton 1 ^. 



III. THEORETICAL SIMULATIONS AND 
PREDICTIONS 

We first clarify that with no phonon coupling included, 
then our Green function technique yields identical nor- 
malized spectra to a master equation approach in the 
low power limit. Specifically, we use the equations in 
Ref.— , with a low-power incoherent exciton pump, and 
we compute the total spectra from both cavity emission 
and QD emission in the presence of pure dephasing; with 
the Green function approach, we define Tzpl = r rac i+r'. 
These spectral forms are found to be identical, as previ- 
ously discussed^, so the general phenomenon of cavity 
feeding is not unique to pure dephasing processes. The 




lu — uj x (meV) 



4-2 2 

lo — lo x (meV) 



FIG. 1: (Color online) Phonon self energies for InAs QDs, for 
two different temperatures (T = 4, 40 K), where the red-dark 
curves represent the imaginary contribution and the grey-light 
curves represent the real contribution. The broadening pa- 
rameters are r ra d = 2 fieV and T' x = 75 /ieV. As discussed in 
the text, lo x is considered to already include a polaron shift. 
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FIG. 2: (Color online) (a) Cavity-emitted spectra for various 
QD-cavity detunings, obtained for a temperature of T = 4K. 
The red-dark and grey-light curves display the Green func- 
tion solution with and without LA-phonon interactions. The 
corresponding imaginary part of the polarizability with (red- 
dark curve) and without (grey-light curve) phonon coupling 
is shown in (b). The lower frequency phonon bath is sup- 
pressed in part due to the low temperatures and in part due 
to the relatively large ZPL. The polaron shift, A ~ 30 /ieV, 
is not included as this just adds in a fixed resonance shift 
for all detunings and temperatures. The parameters are 
r ra d = 2/ieV, u x = 830 meV, T' x = 75 fieV, F c = 100 /xeV, 
and g = 0.13meV (see Ref4£). 



only difference is the overall magnitude, which is due to 
the different initial conditions, though these can be made 
equivalent as well. Any exciton broadening (radiative or 
non-radiative) will feed the cavity mode in the cavity- 
mode emission, and scale with g 2 . However, the distinc- 
tion of pure dephasing plays a much more important role 
for computing higher-order quantum correlations^. 

In Fig. [21b), we show the imaginary part of the QD 
polarizability (which is proportional to the absorption) 
with and without coupling to phonons via the IBM; this 
calculation was obtained at a temperature of 4K, and 
we clearly obtain the familiar spectral form of the IBM 
spectral lineshapei^. The phonon-induced lineshape is 
subsequently used to obtain the modified spectra with 
phonon coupling. In Fig. [2{a) we plot the Green function 
cavity spectra solution with (red-dark curve) and with- 
out (grey- light curve) phonons, where the QD and cavity 
parameters are given in the figure caption. At zero de- 
tuning, we recover an asymmetric Rabi splitting in agree- 
ment with Refsi 13 ' 16 ; in addition, for small detunings, a 
positive detuning gives a cavity suppression, while a neg- 
ative detuning gives a cavity enhancement. These are 
due to the significant frequency shifts from the real part 
of the phonon self-energy, whose effects would be absent 
in an effective rate approximation. This suppression, fol- 
lowed by an enhancement, would exacerbate the effects of 
cavity coupling for positive detunings, and depend upon 
the spectral form of the phonon bath. The off-resonant 
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FIG. 3: (Color online) As in Fig. but with T = 40 K and 
F' = 150 ^V. 



cavity peak primarily stems from S c and scales with g 2 . 

It is also worth stressing that the asymmetric vacuum 
Rabi doublet is an effect that cannot be predicted with 
a standard master equation approach, nor by using an 
effective cavity feeding ratei^. As pointed out some time 
ago by Carmichael and Walls^i, the usual decomposition 
of the system Hamiltonian to include only the nonintcr- 
acting QD and cavity parts does not satisfy the detailed 
balance conditions!. This is due to the fact that the 
cavity and the QD systems are internally coupled, and 
modified master equation must be derived to account for 
this internal coupling. In fact, the polaron transforma- 
tion used by Wilson- Rac and Imamogliiii adopts such an 
approach, and so their transformed system Hamiltonian 
leads to the correct form of the density operator - since 
it preserves detailed balance conditions^. Physically the 
asymmetric vacuum Rabi splitting^, which has recently 
been measured^, comes form the fact that the two cou- 
pled Rabi peaks sample different parts of the asymmetric 
phonon bath. 

Since we include phonons to all orders, it is also pos- 
sible to study elevated temperatures, which is highly 
desired from a device viewpoint. In Fig. [3] we repeat 
the same calculation as before, but at T = 40 K (cho- 
sen primarily to connect to the experiments below), and 
the ZPL has been broadened (doubled) to be consistent 
with experiment^. Again, we see the same phonon cou- 
pling trends as before, but now with a significantly larger 
phonon feeding to the cavity mode. Also, the effect of 
the cavity coupling parameters is significantly reduced 
by about 50%, resulting in a vanishing Rabi splitting on 
resonance when phonons are included. All of these pre- 
dicted trends are consistent with the recent experiments 
of Ota et al^. 
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IV. COMPARISON WITH RECENT 
EXPERIMENTAL DATA ON A SINGLE QD- 
PHOTONIC CRYSTAL CAVITY SYSTEM 

We now apply our theory to help explain the recent 
model discrepancies that were employed for simulating 
experimental data. Earlier this year, Dalacu et al£^ 
reported measurements on a single QD-cavity system 
at low power and, using a master equation theory of 
Refs i i , and clearly demonstrated that an extra cavity 
pump term had to be included by hand to explain their 
data. The cavity pump contribution was found to be 
cavity-QD dctuning-dcpcndcnt. The origin of this cav- 
ity pump was unknown and cited to be somewhat mys- 
terious for their excitation conditions and sample, but 
similar couplings have been shown already to be due 
to electron-phonon couplin g 14 i 16 i 18 . The site-controlled 
QD allows one to suppress other extraneous cavity feed- 
ing mechanism o 21 ' 27 , and for detunings greater than 
5 mcV, they observe no cavity-mode emission at low tem- 
peratures; this is expected for a single exciton-cavity 
coupling 2 ^, since the cavity-mode emission diminishes as 
a function of detuning. 

The single QD-cavity system is realized by nucleat- 
ing one InAs QD at the apex of a InP pyramid (see 
Fig. 2|) grown using sclcctive-area epitaxy^. Dot for- 
mation on these InP pyramidal nanotcmplatcs proceeds 
via the Stranski-Krastanow growth mode^ similar to 
growth on planar substrates although some subtleties of 
the strain distribution will differ due to the proximity 
of the {110} planes that make up the sidewalks of the 
InP pyramid. Once planarized, however, any signature 
of the InP pyramid vanishes, leaving a coherent InAs 
dot in a uniform InP matrix. One distinction between 
planar and site-controlled dots pertains to the wetting 
layer. Although the presence of a wetting layer is as- 




FIG. 4: A single InAs quantum dot nucleated at the apex 
of a InP pyramid. The pyramid is grown using selective- 
area epitaxy on a electron-beam patterned Si02-coated InP 
substrate. The scale bar is 210 nm. 




-2 2 4 -2 2 4 



lu — u x (meV) u — u x (meV) 

FIG. 5: (Color online) (a) Recently published experimen- 
tal data (Dalacu et ai.— ), taken at T = 4K, normalized 
to show the effect of detuning on the exciton and cavity 
mode, (b) Theoretical simulations with (red-dark) and with- 
out (grey-light) phonon interactions. The parameters are sim- 
ilar to those used in Fig. (1), except that V = 150 /ieV and 
we have convolved with a Lorentzian function with FWHM 
r spC c = 250 ^teV to account for the spectral resolution in the 
experiment^. The detunings are not exactly fit to experi- 
ment, but rather are chosen to cover a similar range to the 
experiments with equal frequency spacing. Both sets of cal- 
culations are normalized to their peak value for clarity. We 
also include the radiation- mode emission, with F c = 2F r . 



sumcd, its lateral extent is limited to the apex of the 
pyramid and is thus only slightly larger than the QD. 
The site-controlled dots are therefore expected to have 
a modified shape compared to planar dots, and the ab- 
sence of an infinite 2D wetting layer. Although these 
differences to do not manifest in the electronic structure 
of the dots^i, the absence of a wetting layer and associ- 
ated continuum of states may have implications related 
to non-resonant dot-cavity coupling^ 7 -. Clearly, reducing 
unknown and unwanted excitation mechanisms should be 
avoided from a practical perspective, which makes dots 
without a wetting layer advantageous, e.g., for creating 
cavity-assisted single photons on demand^ 2 -. 

In Fig. 03a) we show the data of Dalacu et al^, and 
in Fig. [SJb) we show our analytical model with (red-dark 
curve) and without (grey-light curve) phonon coupling. 
Evidently, the analytic phonon model well reproduces the 
experiment data without the need to artificially add in 
a cavity pump term, and the phonon coupling plays a 
role on the cavity feeding mechanism. This general con- 
clusion is consistent with earlier work a 15 ' 16 i 18 ' 30 , though 
some of the qualitative trends are quite different (e.g., a 
cavity suppression followed by an enhancement); more- 
over, since our final spectral forms are analytic, one can 
sweep a wide range of parameters in a straightforward 
and simple way. Moreover, our formulas can be applied 
to fit spectra over a wide range of temperatures (up to 
room temperature), where one must go beyond a second- 
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FIG. 6: (Color online) (a) Measured spectra at various tem- 
peratures, taken at T = 10 K, T = 20 K, and T = 40 K. (b) 
Theoretical simulations with (red-solid curve) and without 
(grey-light curve) phonon interactions. The parameters are 
similar to those used in Fig. 5, except that V = 225 — 300 fieV 
(increasing linearly with temperature) and we have convolved 
with a Lorentzian function with FWHM T Bpcc = 400 ^eV— . 



order approximation for the phonon baths. 

Finally, we fit the emission spectra for elevated 
temperatures, using data that has not previously been 
shown. The sample and experiment are similar to 
those described by Dalacu et al£&. In Fig. [5] we show 
the experimental data (a) and simulations (b) for a 
temperature-controller cavity detuning, for temperatures 
of 10 K, 20 K, and 40 K. Once more, we see a good trend 
with the experiments, and the effect of the phonons is 
to enhance the oscillator strength of the cavity mode 
relative to the exciton mode, which becomes stronger as 
a function of temperature. 



We have presented a Green function theory to describe 
photon emission in arbitrary QD-cavity systems, without 
recourse to either the weak coupling regime or a per- 
turbative approximation for the phonons. We have ex- 
ploited this approach to model the strong coupling regime 
between a single exciton and a photonic crystal cavity, 
and, using a very simple initial condition of an inverted 
electron-hole pair, obtained the emission spectra for var- 
ious dot-cavity detunings. Phonon related effects such as 
cavity-mode suppression and enhanced cavity feeding are 
demonstrated, showing the need to include both real and 
imaginary phonon self-energy terms with internal cou- 
pling. The model was then applied to help explain ex- 
perimental data as a function of cavity-exciton detuning 
and as a function of temperature; good agreement has 
been found without artificially adding in a cavity-mode 
pump — as was done previously^. The phonon-induced 
cavity coupling thus mimics an incoherent cavity pump 
term, which naturally would be detuning dependent be- 
cause of the asymmetric phonon bath. Our results and 
general predictions are also consistent with the recent 
data and polaron master equation-^ simulations of Ota 
et alH. 
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